Scaling and energy transfer in rotating turbulence 
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The inertial-range properties of quasi-stationary hydrodynamic turbulence under solid-body rota- 
tion are studied via high-resolution direct numerical simulations. For strong rotation the nonlinear 
energy cascade exhibits depletion and a pronounced anisotropy with the energy flux proceeding 
mainly perpendicularly to the rotation axis. This corresponds to a transition towards a quasi-two- 
dimensional flow similar to a linear Taylor-Proudman state. In contrast to the energy spectrum 
along the rotation axis which does not scale self-similarly, the perpendicular spectrum displays an 
inertial range with fc]| 2 -behavior. A new phenomenology gives a rationale for the observations. The 
scaling exponents £ p of structure functions up to order p = 8 measured perpendicular to the rotation 
axis indicate reduced intermittency with increasing rotation rate. The proposed phenomenology is 
consistent with the inferred asymptotic non-intermittent behavior Q p = p/2. 
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The inherent properties of turbulence in a rotating reference frame are important for, e.g., the dynamics of atmo- 
sphere and oceans, liquid planetary cores, and engineering problems. The nonlinear spectral transfer of energy by the 
direct turbulent cascade and the associated energy spectrum are particularly interesting due to their direct connection 
to the dynamical processes governing rotating turbulence. Most of the available experimental data [3, 0, & BL HI yields 
no conclusive information on the expected self-similar scaling of the energy spectrum in the inertial range of scales 
and its dependence on the rotation frequency fl. Although recent experiments @, 0, Q have shed some light on these 
issues, the scaling of two-point statistics in rota ting turbulence remains a controversial topic. 

Direct numerical simulations [t| E3, EH, E2, EH, Eil 15 , 3| and large-eddy simulations, see e.g. 17, 3, 19| . have been 



carried out only at low and moderate Reynolds numbers precluding clear scaling observations. Nevertheless, most 
of the cited works agree in that the nonlinear spectral transfer of energy to smaller scales diminishes with growing 
fi, accompanied by a transition of the flow towards a quasi-two-dimensional state perpendicular to the fixed rotation 
axis, f2. The transition manifests itself in an increasing ratio of integral length scales parallel and perpendicular to 
f2 = fie 3 , defined as = J^°° d£(vi(r)vi(r + £ej)) / (vf(r)}, representing the largest possible distance between 
two points in the simulation volume and t denoting the respective space increment. 

This Letter presents high- resolution direct numerical simulations of incompressible rotating homogeneous turbulence 
driven at largest scales and proposes a phenomenology of the energy cascade which suggests a physical explanation for 
the observed attenuation of nonlinear spectral transfer under the influence of rotation. In addition, the model gives 
a rationale for the observed trend towards two-dimensionality in rapidly rotating turbulence which is corroborated 
by the simulations. The scaling of two-point structure functions perpendicular to $7 indicates a decreasing level of 
intermittency with growing f2. 

Incompressible hydrodynamic turbulence under solid-body rotation is usually described by the Navier-Stokes equa- 
tions including the Coriolis force [2fj| given here in dimensionless form with the vorticity, uj — V x v, and the 
non-dimensional kinematic viscosity fi and rotation frequency = |f2|, 

dt<jJ — V x (v x uj + 2v x fl) + fiAu) , (1) 
V-i> = 0. (2) 

Equations ((T|) and © are numerically integrated using an explicit trapezoidal leapfrog scheme plj with the diffusive 
term included by an integrating factor, see e.g. [22l |. while the remaining right-hand-side of Eq. ([1]) is determined 
pseudospectrally. The simulation volume extends 2-7r in each dimension with triply periodic boundary conditions and 
a resolution of 512 3 collocation points. Aliasing errors are treated by spherical mode truncation [23^ . Quasi-stationary 
turbulence is generated by a forcing which freezes all modes in a sphere of radius kf = 2. 

The initial state of the forced simulations is taken from non-rotating turbulence which has been freely decaying 
for about one large-eddy-turnover time, the period needed to reach the maximum of dissipation when starting with 
a smooth velocity field. This initial velocity field is characterized by an energy spectrum ~ exp(— k 2 /^), ko = 4, 
and random phases. Subsequently all Fourier modes with k < 2 are frozen. These modes sustain a gentle driving of 
the flow by nonlinear interaction with the fluctuating part of the system. As soon as total energy, E — J v dVv 2 /2, 
and dissipation, e = — /i JydVcu 2 , are statistically stationary with E ~ 1 and e ~ 0.25 both mildly fluctuating, fl is 
set to a finite value, 5 (system I) and 50 (system II). 
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The dimensionless kinematic viscosity /i is 4 x 10~ 4 in both cases. The characteristic length Lq and velocity Vq, 
necessary for the calculation of macroscopic Rossby number, Ro = Vo/(2QL ), and Reynolds number, Re = L Vo/n, 
can only be determined a posteriori in homogeneous turbulence. Both quantities are estimated on dimensional grounds 
using E, e, and 17 as L ~ Ej (lie) 1 / 2 and Vq ~ -E 1 / 2 . Hence the Rossby and Reynolds number given in this paper are 
defined as Ro = y/e/ (4ft E) and Re = yl E 3 / (f2e) / /x, respectively. Note that another common estimate of the Rossby 
number is Ro* = 4Ro 2 . 

After the sudden onset of rotation, E displays a sharp drop of about 20% (case I) and 13% (case II) with a 
subsequent remount that levels off near the previous state. The dissipation rate e follows the general dynamics of 
energy, but does not increase again with e ~ 0.05 in both simulations. The observations can be understood by the 
rotation-induced depletion of the spectral energy transfer which causes a transient until forcing and cascade have 
reached a new equilibrium. The observed behavior does not differ qualitatively if the rotation is ramped up (as has 
been checked by test computations). 

The following 15 (I) and 10 (II) large-eddy turnover times of statistically stationary rotating turbulence are charac- 
terized by Ro ~ 5.3 x 10~ 2 , Re ~ 4000 (I) and Ro ~ 1.3 x 10~ 2 , Re ~ 2340 (II). Perpendicular one-dimensional energy 
spectra, Ek ± = J dk\\ J dk'\v k \ 2 /2 with k± _L ft, || fl, and k' perpendicular to k± and k\\, are shown in Fig. [TJ The 
spectrum of system I displays a scaling range for 4 < kj_ < 20 while in simulation II the dissipation region is starting 
at smaller k (see below). In addition, the spectrum of simulation II exhibits a bump around kf where the forcing 
region descends into the freely evolving range of scales. This effect which is caused by the simple forcing scheme does 
not seem to significantly perturb the flow beyond k ss 5. It is therefore tolerable at the chosen numerical resolution. 

Although the inertial range in simulation II is shorter than in simulation I the perpendicular one-dimensional 
energy spectra in both cases exhibit scaling, Ek ± ~ kj 2 , in agreement with direct numerical simulations at moderate 
Reynolds-number 12] and fl3f ( for k > kf), and shell- model calculations 24 1. A formal analysis of the energy flux in 
helical mode decomposition [25[ leads to the same result as well as dimensional analysis of the energy flux terms which 
occur in quasi-normal closure theories 26|, 27 1 when assuming t* ~ tq ~ for the relaxation timescale of nonlinear 
interactions, r*. The observed spectra are also in accord with one of several isotropic scalings proposed in weak- 
turbulence theory [H, [2^|. However, this approach is only valid in the asymptotic limit t* 3> tq implying k -C kn (see 
below) which requires very strong rotation or an extremely broad inertial range and therefore is beyond the scope of 
the simulations considered here. In [l3| fc _3 -scaling is found at small wavenumbers in turbulence with medium-scale 
forcing. The observation is explained dimensionally by the missing explicit e-dependence of the spectrum at large 
scales (also cf. [lo| using a method based on reduced sets of nonlinear interactions). Recent experimental results Q 
suggest an energy scaling exponent w —2.5 for case II with the micro Rossby number of Ro u = (w|) 1 / 2 /(217) ~ 0.08 and 
an exponent of w —1.7 for case I with Ro w ~ 0.7. This disagreement is probably due to the different way of turbulence 
generation. While in the simulations there is a continuous large-scale forcing, turbulence in the experiment is exited 
initially and then subject to decay under rotation. In addition the experimental Rossby numbers are significantly 
larger than in our computations. The configuration is, therefore, not directly comparable to the simulations described 
here. This is also true for the experiment reported in 0] although the same scaling ~ k~ 2 is observed there (however 
in an inverse energy cascade). It should be noted that the fcii -dependence of the perpendicular energy s pec trum 
(not shown) in our simulations confirms the expected concentration of energy around the plane fcii = [3lL l32j |. 
Furthermore, the perpendicular energy spectra at specific fixed k\\ do not show clear spectral scaling. This is only 
seen in the sum Ek ± ■ 

The energy spectra taken parallel to 1~2 do not exhibit distinct scaling ranges. This is in accord with 
the strong rotation- induced decrease of the axis-parallel nonlinear energy flux (cf. Fig. [2|), TL = 

J^" dfcy J dk 1 f dk± (iu;* • (fe x [v X w] fe ) + c.c^j /k 2 , with [•] denoting Fourier transformation and T)- ± defined anal- 
ogously. Note that the Coriolis force has no direct effect on the kinetic energy since it is oriented perpendicularly to 
v. It does, however, modify the nonlinear interactions leading to depletion and anisotropy of the energy cascade, cf. 
also [jiH [3JI, which is apparent when regarding Fig. [2] All normalized transfer functions for Cl = 0, 5, and 50 are 
negative throughout indicating a direct energy cascade towards small-scales. They show a damping of the energy flux 
at all scales with increasing f2. The cascade depletion is much stronger in than in Tk ± . Consequently, the cascade 
becomes highly anisotropic in the case of strong rotation which naturally leads to a dynamical two-dimensionalization 
of the flow. The transition towards 2D also manifests itself in the increasing ratio of parallel to perpendicular integral 
length scales with growing f2 where strong rotation, Q = 50, yields Li.3 up to a factor 2.5 larger than L^i and for 
i = 1,2,3. This trend is also seen in visualizations of the velocity field (not shown). The perpendicular longitudinal 
and lateral correlation lengths exhibit moderate systematic differences due to the forcing scheme which is not fully 
isotropic. This, however, does not cloud the main trend of two-dimensionalization in planes perpendicular to f2. 
The obtained results can be understood in the framework of a phenomenology of the energy cascade in rotating 
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turbulence which is set out in the following. While inertial waves present in rotating turbulence are undoubtedly 
an important dynamical process, the proposed picture of the energy cascade is based on the presumed energetic 
dominance of convective motions. 

The wavenumber range strongly influenced by rotation lies above the scale at which the advection nonlinearity 
in Eq. is roughly equal to the Coriolis force, v 2 /£ ~ Slve where vi denotes an isotropically estimated velocity 
fluctuation at scale t ~ fc _1 . Since non-rotating hydrodynamic turbulence exhibits Kolmogorov inertial-range scaling 
(33|, ve ~ (rf) 1 / 3 , one obtains the well-known rotation wavenumber 25|, 34 1. kn = (fi 3 /^) 1 / 2 , below which the energy 



cascade is modified by the Coriolis force acting in planes perpendicular to n. 

The nonlinear energy transfer in those planes can be estimated dimensionally as e ~ u| /r tr with the characteristic 
velocity fluctuation in the axis-perpendicular plane at scale £ ~ kj 1 . Due to the lack of an inertial range in the 
direction parallel to f2, a consequence of the quasi-2D-statc, the approximation E(£) ~ u? suffices for the following 
scaling predictions. In isotropic non-rotating hydrodynamic turbulence fluid particles which belong to a turbulent 
structure at scale £ ~ £ follow trajectories of length £ for a turnover time, r tr ~ tnl ~ £A>£j m the course of the 
energy cascade. When neglecting all nonlinearities fluid particles in rotating flows are displaced from the no-rotation 
trajectories in axis-perpendicular directions due to the Coriolis force and follow circular orbits of radius r ~ v^tq 
which close after tq ~ f2 _1 . The nonlinear terms in Eq. fl} which dominate turbulence dynamics cause a strong 
deformation of the circular orbits. The deformed circles which we regard as abstract entities whose shapes can be left 
unspecified will be referred to as 'arcs'. Typically these do not close in themselves after tq but lead to an effective 
axis-perpendicular displacement ~ r of the fluid particles. It is easy to verify that in the rotation dominated range 
of scales, k < kn, the displacement r is always smaller then £, requiring the fluid particles to perform £/r 'arc steps' 
to cover the distance £ and to complete the cascade trajectory. In fact though, the direction in which a fluid particle 
is deflected by an arc movement is quasi-stochastic. Consequently, analogous to a random walk process (£/V) arc 
steps are necessary to cover the distance £. Therefore, the nonlinear transfer time in the axis perpendicular direction 
is given by 

r tr ~ (^) tq . (3) 

This result has also been obtained by formal analysis of the nonlinear energy flux [25| and in weak-turbulence theory 
[28| . We note in passing that the rotation- modified transfer time r tr is larger than t^l with tjvl ~ T tr a t fe- From 
the previous it is clear that for k <C fca the turbulent energy cascade is highly anisotropic and progresses mainly in the 
direction perpendicular to the rotation axis in accord with the dynamic Taylor-Proudman theorem, see e.g. 1|| 27 1. 



The definition of tq used in this paper differs from the one known from weak turbulence theory which involves the 
ratio fcj_/fc|| (see, e.g., [28jj ) . However, both definitions approach each other when k±/kn for the turbulent fluctuations 
does not depart strongly from unity as is the case here. 

With relation ([3]) the dimensional estimate of the nonlinear energy flux is obtained as e ~ v|/(f2£ 2 ). Assuming 
e = const, throughout the inertial range yields a scaling law for the velocity fluctuations perpendicular to ft, 

v e ~(Qe)^e /2 , (4) 



corresponding to the observed scaling of the energy spectrum ~ k J [251 . 1261 ] . 

The rotation-dominated range of scales is limited from below by the wavenumber min(fcn, fcj) where fcj indicates 
approximate equality of nonlinear and dissipative energy-fluxes, e ~ /^£~ 2 u|. Together with relation ((4|) this yields 



~ /^ _1 (e/17) 1 / 2 ~ kQ-^z. A different approach [25(] leads to the same result. Since kn grows and A;^ diminishes 
with increasing f2, there exists a critical rotation frequency for which the rotation-dominated regime has its largest 
extent, i.e. R = e/(/i0 2 ) = I. Higher rotation rates lead to a reduction of this range since R < 1. Here, for case I: 
k n ~ 50, k n ~ 250, k d ~ 167, R ~ 5 and for case II: k n ~ 1581, k n ~ 79, k d ~ 167, R ~ 0.05 which explains the 
shorter inertial range of E^ x in case II. The Kolmogorov dissipation wavenumber k d ~ (s/fj, 3 ) 1 ^ 4 roughly identifies 
the scale at which dissipation begins to dominate over the nonlinear energy flux. 

At the given spatial resolution scaling exponents of the axis-perpendicular longitudinal velocity structure functions, 
(r) — v(r + £)] • ~ can be determined for both systems up to order 8. However, the small inertial range 

for Q — 50 necessitates use of the extended self-similarity property (ESS) [3^|. The relative exponents Cp/C2 obtained 
via ESS approximately coincide with the £ p since relation ^ (in the non-intermittent limit) and the numerical result 
E k± - kj 2 = k~ iC2+1) suggest ( 2 w 1. The results shown in Fig. [3] are a sign of a gradual transition from the 
intermittent non-rotating case (represented by the She-Leveque formula (36|, £ p = p/9 + 2[1 — (2/3) p / 3 ]) towards a 
strongly rotating flow with weak intermittency. The observed reduction of intcrmittency is in accord with the expected 
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FIG. 1: One-dimensional energy spectra perpendicular to the rotation axis for different rotation frequencies. The spectra are 
time-averaged over 6 large-eddy turnover times and compensated with k\ (solid line: Q. — 5, dotted line: S7 = 50). 
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FIG. 2: Normalized and time-averaged nonlinear energy fluxes (a) parallel and (b) perpendicular to the rotation axis with 
Q, = (dashed line), Q, = 5 (solid line), and Q = 50 (dotted line). 

transition from strong fluid turbulence to weak inertial wave turbulence (see, e.g., [371 )• The observations coincide 
with experimental findings given in [7] although in this experiment the flow is driven at small scales and in contrast 
to our simulations exhibits an inverse cascade of energy. 

In summary, high-resolution direct numerical simulations of incompressible hydrodynamic turbulence driven at 
largest scales under moderate and strong rotation corroborate a proposed phenomenology of rotating turbulence 
which gives a simple rationale for the overall weakening of the energy cascade and the trend towards two-dimensional 
dynamics in rotating turbulent flows. Higher-order structure function scalings show a transition towards a non- 
intermittent state perpendicular to the rotation axis as known from 2D turbulence. 

It is a pleasure to thank Friedrich Busse and Dieter Biskamp for helpful discussions. 
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FIG. 3: Axis-perpendicular structure function scaling exponents ( p of velocity up to order 8 for moderate and strong rotation. 
Solid line: Q = 5, dotted line: f2 = 50, dashed line: She-Leveque intermittency model, dash-dotted line: non-intermittent 
scaling, ( p = p/2. 
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